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ABSTRACT 

The eclipsing binary X-ray source 4U 1700 — 37 has been continually monitored by the BATSE experi- 
ment on the Compton Gamma Ray Observatory since the spring of 1991. Using source measurements at 
times of Earth occultation, we observe an average (uneclipsed) flux of 0.23 crab in the 20-120 keV band. 
The flux is highly variable, with occasional flaring behavior on timescales from hundreds of seconds to 
several hours and intensities as bright as 1 crab. The uneclipsed spectrum is well represented by an opti- 
cally thin thermal bremsstrahlung model with a temperature of 25 keV independent of source intensity 
or orbital phase. An upper limit of 4% on the pulse fraction has been obtained for pulse periods between 
2 and 700 s. Average orbital light curves from almost 1000 days of occultation measurements have been 
constructed. These profiles are used to measure (1) the eclipse semiangle, 0 E = 28?6 + 2?1 in the 20-120 
keV band, and (2) the decrease in orbital period, P/P = -(3.3 ± 0.6) x 10 -6 yr' 1 . Estimates of system 
physical parameters are obtained using Monte Carlo simulations to propagate errors in measured and 
assumed parameters. For the X-ray source mass we find M x = M 0 , and for the mass and radius 

of the optical companion, M 0 = 30! 7 1 Af 0 and R 0 — 18 tj R Q . 

Subject headings: binaries: eclipsing — stars: individual (4U 1700 — 37) — stars: neutron — X-rays: stars 


1. INTRODUCTION 

The high-mass X-ray binary (HMXB) 4U 1700 — 37 was 
first detected with the Uhuru satellite (Jones et al. 1973). 
Uhuru observations revealed that it is an eclipsing binary 
with a period of 3.412 days. Optical observations of the 
companion HD 153919 indicate that it is an 06f star, 
located at a distance of about 1.8 kpc. From these observa- 
tions it was inferred that HD 153919 slightly underfills its 
Roche lobe (Conti 1978). This implies that the X-ray source 
is powered by the wind of the companion star. The mass of 
the X-ray source is consistent with that of a neutron star 
(Hutchings 1974), and the high-energy spectrum is similar 
to that of X-ray pulsars (Dolan et al. 1980). This suggests 
that 4U 1700—37 is a neutron star, despite the fact that no 
X-ray pulsations or X-ray bursts have ever been confirmed. 

The lack of a confirmed pulse period complicates the 
determination of system parameters. This determination is 
further limited by the uncertainty on the companion radial 
velocity semiamplitude K 0 and the duration of X-ray 
eclipse. Here we will reexamine the system parameters 
based on an improved eclipse semiangle measurement. 

The decrease in the orbital period, first observed by 
EX OS AT (Haberl, White, & Kallman 1989), is confirmed 
by our observations. Two possible causes of orbital period 
decrease in this system are those due to wind-driven 
angular momentum loss from the companion star and the 
tidal mechanism. By comparing our result to the expected 
rates of change from these sources, we find that angular 
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momentum loss is most likely the dominant cause of the 
period change. 

The BATSE experiment has obtained the longest contin- 
uous record of observations of 4U 1700—37, sampling its 
flux up to 30 times per day for over 1000 days and recording 
flares at 1.024 s resolution. In this paper we will examine the 
characteristics of the system exhibited in the long-duration 
BATSE database. In § 2 we present the observations. The 
eclipse semiangle is determined by fitting the orbital profile 
of the X-ray emission to a model. Using this result, the 
system physical parameters are calculated in § 3. The 
orbital period change is determined in § 4. These results are 
further discussed in § 5, and our conclusions are sum- 
marized in § 6. 

2. OBSERVATIONS 

The Burst and Transient Source Experiment (BATSE) on 
the Compton Gamma Ray Observatory (CGRO) is an all-sky 
monitor for gamma-ray bursts, persistent sources, and tran- 
sient phenomena. Here we report on the analysis of data 
obtained with eight uncollimated Nal(Tl) large area detec- 
tors (LADs) (effective area ~ 1500 cm 2 per detector in the 
20-50 keV band, sensitive to photons in the 20-1800 keV 
band) which are mounted in an octahedral geometry 
around the CGRO spacecraft (Fishman et al. 1989). Data 
with 1.024 s time resolution is telemetered to the ground in 
four energy channels (DISCLA data type), and with 2.048 s 
resolution in 16 energy channels (CONT data type). 

About 85% of the sky is occulted by the Earth twice every 
spacecraft orbit (~90 minutes). When a source is sufficient- 
ly bright (>300 mcrab), the profile of Earth occultation is 
visible in the raw count rate data. At each Earth 
occultation, count rate measurements are accomplished by 
separately fitting each energy channel of CONT rates to a 
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model. The model consists of a quadratic background 
model plus terms for known bright sources in a 4 minute 
window centered on the limb crossing time of the source 
(Harmon et al. 1992). We call these edge fits. Each detector 
with its axis less than 60° to the source (between one and 
four detectors, but typically two) is fitted separately. The 
source-term profile as a function of time is determined from 
an energy-dependent model of the atmospheric transmis- 
sion along the line of sight to the source. Typically, the 
transmission varies between null and full during a 20 s inter- 
val, but this interval can vary from 10 to 90 s depending on 
the angle of the source to the plane of the spacecraft orbit. 
Photon fluxes can then be obtained by convolving a spec- 
tral model through detector response matrices and fitting to 
the measured count rates. The 3 o sensitivity obtained by 
averaging occultation edge measurements is about 100 
mcrab per day. 

The quadratic background model used in the edge fits is 
inadequate for windows longer than about 300 s. In order 
to represent detector background effectively on longer time- 
scales, a global background model has been developed 
(Rubin et al. 1992). This model contains terms which rep- 
resent backgrounds due to the diffuse cosmic background, 
prompt radiation caused by cosmic-ray interactions in the 
atmosphere and spacecraft, radioactivation of detector 
materials during South Atlantic Anomaly (SAA) passages, 
and bright occulting sources. The occulting source trans- 
mission model used in the edge fits is also used in the back- 
ground model. Below we use background-subtracted 
DISCLA data to obtain light curves and perform timing 
analysis of flaring episodes. Background-subtracted data 
have also been used to search for coherent pulsations from 
4U 1700-37. 

2.1. Light Curves of Flaring Events 

Flaring behavior from 4U 1700-37 is observed as an 
increase in count rate visible in the unocculted LADs facing 
the source. Changes on timescales from minutes to hours 
are observed in the DISCLA data. Isolated flaring events 
are observed, as well as series of flaring events lasting from 
hours to 2 or more days. That an increase in the back- 
ground count rate is associated with a flare on 4U 1700 — 37 
can be verified by using the relative count rates observed in 
different detectors to determine a sky location for a flaring 
event. Algorithms previously developed for locating 
gamma-ray bursts (Brock et al. 1992) have a statistical 
uncertainty of about 15° for untriggered events with signifi- 
cant emission between 20 and 100 keV. If several bright 
occultation steps are available, these can be used to verify 
the source location to within 0?5. 

Often a series of flares erupts over 0.5-2 days. A good 
example occurred while the source was uneclipsed in its 
binary orbit from TJD 8904.2 to TJD 8906.8. Figure la is a 
light curve of occultation count rates measured in DISCLA 
channel 1 (20-50 keV) data in LAD 3, which had a viewing 
angle of 19°. In this figure, occultation count rates deter- 
mined from the background modeling are shown. The 
source is quite active during this interval, and variations by 
factors of 2-5 on timescales of hours are not uncommon. 

These hourly variations are accompanied by rapid 
changes on shorter timescales. Figures \b-\d show three 
episodes of fast variations, observed in the background- 
subtracted data of DISCLA channel L During this time, the 
background count rate varies between 800 and 1500 counts 


s 1 . In the first episode shown in Figure 16, the count rate 
due to 4U 1700 — 37 remains essentially constant at a level 
of about 1.5 crab for 700 s, except for one spike of less than 
100 s duration. The location of this spike is within 5?0 ± 6?2 
(0.8 <s) of 4U 1700 — 37. In contrast, the second episode, 
shown in Figure 1 c, contains large, incoherent fluctuations 
on timescales above 50 s. In the third episode, in Figure Id, 
the count rate has only relatively small variations for 2500 s 
before increasing sharply and pulsing at least 3 times. A 
fourth pulse is in progress as the source is occulted by the 
Earth. The pulse peaks are separated by 75 s and have been 
located to within 25?5 + 1 6?6 (1.5 <r)of4U 1700 — 37. 

As discussed in § 2.5, no clear evidence of a pulsar period 
for 4U 1700 — 37 has been found. However, in addition to 
the previous episode, another pulsing episode, this one con- 
taining six intense pulses, has been observed. The flare con- 
taining these pulses occurs at the peak of a series of flares 
lasting about 0.5 days. Background-subtracted DISCLA 
count rates from LAD 2 recorded during this episode are 
shown in Figure 2. The first pulse begins at about 81,950 s, 
and the pulses continue for about 550 s. The time between 
successive pulse peaks is about 100 s for the first five pulses 
and somewhat longer, about 150 s, before the last pulse. No 
pulses are visible in the 300 s before the source sets, 
although the flare is still clearly in progress. The location of 
the third pulse, which is the brightest, is within 20?3 + 10?9 
(1.9 a) of 4U 1700—37. The setting occultation step occurs 
where a setting step from 4U 1 700 — 37 is expected. 

Few nearby sources which could potentially be confused 
with 4U 1700 — 37 exhibit the flaring behavior described 
above. Scorpius X-l flares are visible at lower energies (in 
the BATSE spectroscopy detectors), but rarely above 20 
keV (McNamara, Harmon, & Harrison 1995). The loca- 
tions determined for the flaring events are inconsistent with 
the location of Scorpius X-l at the 3 o level or higher. 

2.2. Long-Term Light Curve 

The long-term light curve from the occultation analysis is 
shown in Figure 3. Each data point is the average of the 
20-120 keV flux during binary phases 0.094-0.906 of one 
binary orbit. The horizontal dashed line is at the mean flux 
of 0.23 crab = 0.069 photons cm -2 s' 1 (hereafter y cm~ 2 
s _l ), corresponding to a luminosity of 1.5 x 10 36 ergs s"L 
The uncertainty of the mean is 0.02 crab. There are no clear 
long-term trends; however, the strong variability of the 
source is apparent. The rms scatter of the flux measure- 
ments is 0.24 crab. For comparison, BATSE occultation 
measurements of the Crab flux vary typically by 5%, due 
mainly to systematic effects. 

2.3. Spectra 

The spectrum of 4U 1700 — 37 has been previously mea- 
sured at low energies, in the 0.5-4 keV range with Einstein , 
the 2-60 keV range with HEAO 1 A-2 (White, Kallman, & 
Swank 1983), and the 2-40 keV range with EXOSAT 
(Haberl et al. 1989). Higher energy measurements, above 20 
keV, can be fitted with an optically thin thermal brems- 
strahlung (OTTB) model. This spectrum was observed on a 
balloon flight with kT ~ 28 keV (Pietsch et al. 1980), with 
OSO 8 with kT ~ 27 keV (Dolan et al. 1980), and with 
SIGMA with kT ^ 45 keV (Sitdikov et al. 1993). In each of 
these observations the shape of the source spectrum 
remained constant, independent of intensity or binary orbit 
phase. Recent observations with the Mir Kvant observatory 
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Fig. 1. — Flaring series observed during TJD 8904-8906. Data from the 20-50 keV band in LAD 3 are displayed, (a) Occultation count rates measured 
with the background model. ( b-d) Background-subtracted count rates of flaring episodes, plotted at 8.2 s resolution. Data shown are count rates which may 
be attributed to 4U 1700 — 37. The statistical error on each data point is in the range 10-15 s. Gaps in the data are due to gaps in the telemetry, (b) Source 
rises at 927 s on the horizontal axis, (c) Source rises at 160 s and sets at 3480 s on the horizontal axis, (d) Source sets at 3530 s on the horizontal axis. Note the 
three pulses separated by ~ 75 s just before set. 
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Fig. 2. — Flaring event on TJD 8552. Residual count rates from the 
20-50 keV band in LAD 2 are shown. The first pulse begins at about 
81,950 s, and the pulses continue for 550 s. The sharp decline at 82,785 s is a 
setting edge of 4U 1700 — 37. 
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Fig. 3. — Long-term light curve. Each data point represents the average 
20-120 keV uneclipsed flux in one binary orbit from the occultation 
analysis. Flux is averaged over binary phases 0.094-0.906. The dashed 
horizontal line is at the mean flux of 0.069y cm “ 2 s” 1 . Points without error 
bars at zero flux represent binary orbits in which no useful measurements 
were made. 
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in the 4-100 keV band have suggested thermal Comp- 
tonization of blackbody photons by a hot plasma as the 
emission mechanism (Maisack et al. 1994). The parameters 
of the plasma are kT ^ 30 keV and optical depth t ~ 2. 

BATSE observations are consistent with these earlier 
results. We observe an OTTB model spectrum with a tem- 
perature of about 25 keV during all uneclipsed phases and 
source intensities. Figure 4 shows the uneclipsed spectrum 
and model averaged over a 2 week spacecraft pointing 
period. In this and other spectra, count rates from six chan- 
nels spanning 20-120 keV were input to the fitting pro- 
cedure. To make this spectrum, 166 detector steps between 
TJD 8450.0 and 8463.0, falling in binary phases 0.109 
0.875, were used. Two detectors were fitted simultaneously. 
The resulting amplitude is A — 0.223 + 0.014y (cm 2 s 
keV)~ 1 and the temperature kT = 25 ± 4 keV. The reduced 
X 2 is x 2 / v = 9.6/10. For other subsets of data, the measured 
temperature is generally consistent with kT = 25 keV. The 
corresponding fits y 2 & re also close to unity. 

2.4. Orbital Profile 

Folding flux and spectral measurements from many 
orbits into bins of orbital phase permits construction of an 
average orbital profile. Both Granat/ watch (Sazonov et al. 
1993) in the 8-60 keV range and Prognoz 9 (Kudryavstev & 
Svertilov 1991) in the 10-100 keV range have averaged con- 
tinuous observations over many binary orbits (about 73 
days and 13 months, respectively), but with low spectral 
resolution. Figure 5 shows the 20-120 keV BATSE flux 
folded into 128 phase bins. Both eclipse egress at about 
phase 0.1 and eclipse ingress at about phase 0.9 are clearly 
visible. Variations of flux with phase are due partly to 
flaring behavior of the source. 

The average phase, <f> p , at which 99% of the source flux is 
occulted at eclipse ingress or egress, and the transition 
phase width t p , were estimated by fitting the 20-120 keV 
flux profile to the empirical function (Clark, Minato, & Mi 
1988): 

Fj = F u exp {In (0.01) exp [-(<£/ - <A P )A P ]} . (1) 

where p = i,e for eclipse ingress or egress, and j is the phase 
bin number. F u , the average uneclipsed flux, was determined 
prior to the fit. The first half of the data was included for 
ingress, and the second half for egress. 
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Fig. 4. — Spectrum of 4U 1700-37. Data points represent average flux 
observed during TJD 8450-8463 in binary phases 0.109-0.875. The solid 
line is the OTTB model spectrum fitted to the data. The model tem- 
perature is 25 ± 4 keV. 



Orbital Phase 

Fig. 5. — Binary orbit profile of 4U 1700 — 37 in the 20-120 keV range 
The data were folded into 128 phase bins at the binary orbital period. The 
error bars during eclipse have been uniformly adjusted to reduce system- 
atic errors in the occultation analysis. Error bars outside of eclipse have 
been further (uniformly) adjusted for consistency with the average 
uneclipsed flux. Data before and after phase 0.5 were fitted separately 
according to the model equation in the text. The solid curve shows the 
best-fit model. 

To account for systematic errors in the occultation 
analysis and flaring behavior of the source, the error bars 
were uniformly adjusted before fitting, in the following 
manner. First, data points during X-ray eclipse were fitted 
to a constant flux value of zero. The error on each of the 
data points was multiplied by the square root of the reduced 
X 2 of this fit. This factor adjusts for systematic error in the 
occultation measurements. Second, the uneclipsed data 
points, after their errors were multiplied by that factor, were 
fitted to a constant. The error on each uneclipsed data point 
was multiplied by the square root of the x 2 of this fit to 
adjust for the flaring. Both systematic and flaring correction 
factors were typically between 1.0 and 1.5. 

The model resulting from the fitting procedure is shown 
in Figure 5. The sizes of the error bars include the adjust- 
ments discussed above. The values of the fitted parameters 
in equation (1) are listed in Table 1. In constructing the 
orbital profile of Figure 5, we have used the ephemeris 
derived from BATSE data in § 4, including the orbital 
period derivative. The sensitivity of the profile parameters 
to the ephemeris has been checked by repeating the fitting 
procedure for periods ranging from 3.411401 to 3.411761 
days in steps of 0.00003 days. This did not lead to signifi- 
cantly smaller transition widths than the values listed in 
Table 1. 

The duration of X-ray eclipse can be calculated from the 
measured phases of ingress and egress. This duration is 
0.159 + 0.012 in phase, which is equivalent to 0.542 ± 0.040 


TABLE 1 

Orbital Profile Fit Parameters 
in the 20-120 keV Band 8 


Parameter 

Value 

*. 

0.0883 T 0.0075 



0.0122 + 0.0038 

*, 

0.9295 -1- 0.0091 

T,- 

0.0271 + 0.0051 

xi/v 

X, ? /v 

62/64 

46/64 

a The parameters are in units of 
the orbital phase. 
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?n^9n h J S v m u Plle i 3n eclipse sernian 8 | e of 28?6 ± 2?1 in the 
20 120 keV band. Evidence for a trend in ecliDse duration 

versus energy has been searched for by dividing the 20-120 
eV range into six bands and separately applying the above 
procedure ,n each band. There is no statistically ^significant 
rend with energy. The OTTB model temperSe “so 
shows no dependence on orbital phase, either near eclipse 
or at any time during the orbit. P 

The BATSE measurement represents more than a factor 
hLh p ed 0n 10 the uncertaint y over the two previous 

21 ° ± n 10 ° by 050 5 (Dola " et 
lov lQOn w 1 b y Prognoz 9 (Kudryavstev & Sverti- 
, d l' lS consistent WIt h the measurement of 25° -f 6° 
9781 Hammerschlag-Hensberge, & z„fd e ™i* 

1978) made by comparing the mideclipse time obtained 
from an optical light curve with the X-ray eclipse egress 

&S5ST by Copermcus (B “’ ETS 

to Uf rf?f th ° f the ing [ ess and egress transitions are related 
o the distance over which source flux is occulted by dense 
matter. However, it should be noted that equation (1) was 

fo^which a£ t0 ! nd u Vid u al binary eclipse transitions 
tor which a single scale height applies. Folding data over 

™ any transdions involves the superposition of differing 

scale heights, and, in general, the measured widths therefore 

wHo": S'l F PhySiC s al r er P:e‘ a tion. With this proviso! 
hrJT !l gUre 5 that ecllpse egress appears to be 

EXoZr 'mTT Tb * differe "“ which w as observed by 
. v (Haberl et al. 1989), is significant in the 20-120 

»-30 £fvt'„ h 5 V r " ,Cl ' The differe "“ 1S Sreateriu !he 
Thi !! «• V b d ’ but 1S not S1 8 ni ficant at higher energies 

paniou w'ntT ^ CaUSed b> an i" the com- 

2.5. Pulse Period Search 

Three previous experiments have claimed observations of 
-ray pulse periods: Copernicus (Branduardi et al 1978) at 
24 mmutes 5 (Matilsky, La Sala, & Jessen ?978) at 97 
minutes, and Tenma (Murakami et al. 1984) at 67 4 s^ Sub 
sequent observations by Copernicus (Branduardi et al. 1979) 

that ?he sTv ? r nh ,r f ° rigina ' o^ervation. It was shown 
the ^^ Si . b V n t,on Probably arose from the effect of 
he 100 minute satellite period and telemetry dropouls 
^ammerschlag-Hensberge, Henrichs, & Shaham 1979) 
flare 67 4 i S PeriOC ! WaS observed with Tenma during one 

shoSnTn S 2 1 S Th ar *° th f f amples fro m BATSE data 
shown in § 2.1 These periods have not been observed at any 

other time. It is therefore unlikely that they represent a 

fud? Tr SUrement ° f the pulse period ' The P uls sd ampli- 
]9 o 3 ™l r f 0m one *j E A° I A-2 observation (White eud. 

SirMA^ 0 f ° r Penods , between 160 ms and 480 s. Recently 

~ a TV et aL 1992) used 1 da y of data to set a 
hmit on the pulse fraction of about 0.3% on pulse oeriods in 

the neighborhood of each of the three previous claims 
Two different techniques to search for coherent pulsa- 
Oons m the BATSE data have been attempted. The first 

S i.7? T CS th i C Ph ° t0n flux measur ements derived from 
were OCca tatlon ed ges. The edge measurement times 
were corrected to the solar system barycenter. A Lomb 

fmm nnf r m T* on 4092 ste P measurements 

rom uneclipsed portions of the light curve This is neces 

sary because these fluxes are sampled unevenly in Le No 

evi ence for pulsations due to the source were found. This 
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implies an upper limit on the pulse fraction of about 20% 
for periods between 300 s and 1 00 days 

Second, the DISCLA data are routinely searched for 
oherent pulsations by fast Fourier transform (FFT) tech- 

c £7slT<,° “f rSE da, i y m ° ni,ori "* The LeuLofa 
clear signal detection in the direction of 4U 1700-37 

implies an upper limit on the pulsed flux of about 4% of the 
average source flux determined from the occasion mea- 
surements for periods between 2 and 300 s. This technique 
becomes less sensitive for periods greater than 300 q s 
2™%*: subtraction of an orbttal background model from 

300 and lh< V*" sitivi '>' for P ul * Periods between 
300 and 700 s. A search of two sets of 3 days of data each of 

which contains bright outbursts, was attempted. No signal 
was observed during these intervals or during 1 mon h o 
routine mouttoring. The upper limit on the p„™d flux 
remains at about 4% for a period of 700 s. 

2.6. Source Variability 

hiih r° rder t0 VieW the source var ' a bility qualitatively at 
high frequencies, we have constructed the average power 
spectrum of background-subtracted DISCLA channel 1 
data (20-50 keV) during TJD 8904.2-8906.8 The TSl im 

which a , ZT ° f 102 1 024 X 512 s M*** iutervalsduring 
S seH« „7n^ Wa n u "°“ u |! ed These intervals contain 
“T, shown Ftgure I. As can be seen ,n 

igure 6, the power begins to rise below 40 mHz and con 
tmues to increase below 10 mHz. There are no persistent 
quasi-penodic oscillation (QPO) peaks present Araverage 

n00 r 37 PeCtmm ° btained durin g intervals when 4U 
1700 - 37 was occulted by the Earth did not show anv 

excess, remaining consistent with zero between 5 and 500 

The power spectrum of the pulsing episode of Figure 2 is 

PU 'T d oSn are 

associated with the enhanced power around 10 mHz Ten- 
tatively, we may call this event a transient QPO with fre 
quency ~10 mHz and width ~5 mHz. The three pulses 
separated by 75 s in the flare of Figure Id mav he 7n 

example of the same type of transient phenomenon at a 
somewhat different frequency. ^nomenon at a 


10.00 
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100 

Frequency (mHz) 

ban?obs7^ - 20-50 keV 

of the rms intensity per frequency intervatarml/^Hz 6 '] 8 Jhere^he 
average intensity .s the average occultation count rate observed durint he 
interval. Pos.t.ve detections are indicated by the histoera ^ ZnA V 8 

•024 x 512 s intervals. Frequencies above 50 mHz have been rebinned 


RUBIN ET AL. 


Vol. 459 


264 



An estimate of the brightness distribution of the source 
can be obtained by making a histogram of the number of 
occultalion edge measurements made (unctran of flu ^ 
The measured distribution can be mode y 
an Experimental resolution function w.th a model for the 
intrinsic flux distribution and then fitting to the measured 

diS Se U flux n distribution of the occultation edge measure- 

account for measured flux values below zero, is apparent 
During intervals when the source is in binary eclipse, th 

ex^Erinsic flu, probabiht, d.trtbuhon ts equa .o . 

delta function at zero flux. Using a Gaussian for the expe 

mental resolution function and fitting to those datayieldsa 
menra .yv, no2v cm ~ 2 s 1 and a standard 

Gaussian mean of 0.004 ± O.OUZj cm s _ 

dC AttemptingEcf fU ^uneET.psed flux distribution to the 
convolution of a Gaussian resolution function with a delta 
function at an intrinsic flux determined by the fit leads to an 



sTThis isEkarly an^unaEcqD table result. Instead, an expo- 
nLSiSS cal be used for .he in.rinsie flu, d,s.„bu„on. 

P{F) = Ae xF , (2) 

where F is the intrinsic flux and A is a no rmahzation con- 
Ttant This yields a standard deviation of 0.535 ± 0_0008 
for the Gaussian resolution function and a va ue of a = 

>5 00 + 0.27 for the exponential constant with X / 

177 5/88 The standard deviation is consistent with 
™ ired for fte e/hpsed dis.nbu.ion, and the X 1 * "'$■ 
buT^ot extraordinarily so. Th.s model ,s shown as the 

^mvX'of^’oSned when fining all of .he uneclipsed 
data to an intrinsic exponential flux model could be high 
because that model is not a good representation of the data 
o^because the resolution function used is not adequate^ 
From Hgure 8 we note that there is an apparent excess of 
flux measurements around the peak. It is thus P oss ‘ b1 ^ 
the data are consistent with an exponential model for fluxe 
above The peak. bu. .ha. .here is an exees above ^he expo- 
nential prediction at low flux values. To test this poss bil v 
the exponential model fit was repeated using only the 
uneclipsed data with fluxes above 0.085y on • 

results in a Gaussian 

«,ive?y'^d“ha‘"?he disuibu.ion of flaring episodes 
ihnve ~01v cm 2 s 1 is consistent with an exponential 
model, but that the source also has a persistent, less intense 

Sta Occultation edge measurements have been used to search 
for eGdenrof correlation of spectra with intensity First 
measured edge count rates during the unechpsed poj 
thp lieht curve were fitted to an OTTB model ana inc 
resu ltfng fluxes ranked from brightest to dimmest. SeconcE 
average fluxes and temperatures fitted by this model were 

te a few^is lower fo, fluxes Jus, above .he mean m 

comparison to more intense measurements, but this 
low significance. 
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From the data presented in this and the previous sub- 
sections, we can establish the following conclusions about 
the flaring behavior: (1) flares do not occur preferentially at 
particular orbital phases; (2) the spectra or temperature of 
flares is independent of intensity; (3) significant variations 
occur on timescales greater than 30 s; and (4) the distribu- 
tion of flare sizes above ~0.1y cm' 2 s _1 is consistent with 
an exponential model. 

3. SYSTEM PARAMETERS 

The physical parameters of at least six binary systems 
containing X-ray pulsars, including Vela X-l and 4U 
1538-52, have been reasonably well determined from 
independent measurements of the orbital period, optical 
radial velocity semiamplitude K a , X-ray eclipse semi- 
angle 9 El and projected semimajor axis of the X-ray source 
a x sin i. Observations of optical ellipsoidal variations were 
used to limit the possible range of the Roche lobe filling 
fraction of the companion (Rappaport & Joss 1983 and 
references therein). For these systems, the effects of uncer- 
tainties in the measured parameters have been propagated 
through formulae relating those parameters to M x , M 0 , R 0 , 
i, and a. This has been accomplished by performing Monte 
Carlo simulations in which independent sets of measured 
parameters are randomly generated according to probabil- 
ity distribution for each of the derived parameters 
(Rappaport & Joss 1983). 

This technique has not previously been applied to 4U 
1700-37 because without a pulse-timing analysis, it has 
not been possible to measure a x sin Without an indepen- 
dent measurement of this or some other parameter, the 
physical parameters of this system are underdetermined. 

Previous analyses (Hutchings 1974; Wolff & Morrison 
1974; Heap & Corcoran 1992) have proceeded by making 
assumptions about the possible values or ranges of the 
orbital inclination i and the O star radius R a . In particular, 
it is sometimes assumed that, because the source is eclips- 
ing, i is close to 90°. However, there exists no independent 
evidence for the validity of this assumption, and values of i 
closer to 60° are just as plausible. There is some information 
available about the inclination of this system from measure- 
ments of the variable optical linear polarization of HD 
153919 (Dolan & Tapia 1984). These measurements yield a 
90% confidence interval of 71° < i < 86°. Previous esti- 
mates or ranges for R 0 have sometimes been derived from 
radiation-driven stellar atmosphere models and have not 
been independently corroborated. It is possible to set 
model-independent lower and upper limits to R 0 from 
measurements of the effective temperature T e and the 
optical luminosity using the relation L 0 = 4?rR£ cr7+ where 
a is the Stefan-Boltzmann constant. Using the extreme 
values implied by the ranges 32,000 < T e ( K) < 40,000 and 
1.2 x 10 39 < L 0 (ergs s“ ') < 4.9 x 10 39 (van Paradijs et al. 
1978) leads to 12 < R 0 < 40 for R a in solar units as a firm 
range for R 0 . 

Some previous analyses (Conti 1978) have used values of 
8 e derived from lower energy experiments. However, 
because the flux from the source is attenuated by a wind, 
there should be a decreasing trend of measured eclipse 
duration with increasing energy. The observations, listed in 
Table 2, are consistent with this expectation. Moreover, 
lower energy measurements, and measurements averaged 
over only one or a few orbital cycles, are susceptible to 
systematic error due to the highly variable nature of the 


TABLE 2 


Eclipse Semiangle and Mideclipse Epoch Measurements 


Experiment 

Energy Range 6 e 

Mideclipse Epoch 
(JD -2,440,000) 

keV 

Degrees 

Vhuru 

2-6 

58 ±4 

1453.14 + 0.01 

Copernicus 

2.8-8.7 

32 + 4 

2609.75 + 0.01 

Copernicus 



2613.146 + 0.01 

Copernicus 


40 + 4 

3002.104 + 0.01 

Copernicus 



3005.5 ± 0.01 

Copernicus 




and Optical 


25 + 6 


OSO 8 

21-84 

21 + 10 


Prognoz 9 

10-50 

26 + 5 


EX OS AT 

2-10 

44.5 ± 0.6 

6161.34 ± 0.0030 

Gramzf/WATCH 

8-20 

37 ± 3 

8723.440 ± 0.031 

BATSE 

20-120 

28.6 ± 2.1 

8651.865 ± 0.031 

BATSE 



9149.925 ± 0.027 


source. The source profile is not smooth, making it diffificult 
to fit, and the flux can fall to the eclipsed level even though 
the source is unocculted, leading to an overestimate of 
eclipse duration. Thus higher energy measurements aver- 
aged over many orbital cycles should give the most accurate 
estimate of 0 E . 

Here we will reexamine the system parameters based on 
the improved eclipse semiangle measurement using BATSE 
data. Monte Carlo simulations will be used to estimate 
parameter probability distributions under two different sets 
of assumptions. The first set will be optimized for measuring 
the X-ray source mass. In the second set, particular values, 
or a range of probable values, of M x will be assumed, and 
the remaining parameters estimated. 

The companion Roche lobe filling fraction Q is defined by 

Ro = OR L , ( 3 ) 

where an approximate formula for the Roche lobe radius is 
(Rappaport & Joss 1983) 

RJa = A - B log q + C(log q) 1 ; (4) 

here a is the semimajor axis of the orbit, and q = M 0 jM x . 
The coefficients A, B, C are defined by 

A = 0.398 - 0.026T 2 + 0.004r 3 , 

B = -0.264 + 0.52T 2 - 0.01 5 r 3 , 

C = -0.023 - 0.005T 2 , 

where T = cj 0 /(o is the ratio of the rotational angular fre- 
quency of the computation star to its orbital angular fre- 
quency. For a corotating system, T = 1. For 4U 1700-37, 

T is in the range 0.5 < T < 1.0 (Conti & Cowley 1975; Hut- 
chings 1974). Since HD 153919 does not overfill its Roche 
lobe, Q < 1. A lower limit for D can be set by noting that 
observed ellipsoidal brightness variations in optical light 
curves imply D > 0.8 (van Paradijs et al. 1978). 

The O star radius is related to the orbital semimajor axis, 
inclination, and eclipse semiangle by 

Ro/ a = ^/cos 2 / + sin 2 i sin 2 W E . (5) 

The companion mass function, in solar units, is determined 
from the measured radial velocity semiamplitude K 0 (in km 
s 2 ), the period P (in days), and the assumption of a circular 
orbit : 

fo= 1.038 x 10' 7 X 3 P, 


( 6 ) 
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where f a is defined by 

M 3 sin 3 i =f 0 (M 0 + M x ) 2 . (7) 

The upper limit on the orbital eccentricity, obtained from 
orbital fits to UV radial velocity curves, is about 0.01 (Heap 
& Corcoran 1992). The semimajor axis of the binary orbit is 
related to the component masses by Kepler’s third law : 

a 3 = 75.19(1 + q)M x P 2 , (8) 

where the masses and binary separation are expressed in 
solar units and the orbital period P is in days. 

Depending on how these equations are combined, it is 
possible to obtain solutions for different combinations of 
input parameters. For example, by combining equations 
(3)— (5), a numerical solution for q can be determined for 
each set of generated values of i, fi, T, 9 E , and K 0 . Equation 
(7) then determines M x . However, this would necessitate 
using an input distribution for i which contains a large 
uncertainty. Values of R 0 determined on output by combin- 
ing equations (5) and (8) could be inconsistent with the 
range derived from optical observations. Applying a con- 
straint on R 0 after the simulation, and retaining only 
parameter sets which satisfy it, would have the effect of 
altering the original empirical probability distributions. In 
general it is desirable to minimize the effect of such 
alterations. 

Another procedure would be to generate values for R 0 , 
Cl, T, d E , and K 0 . Equation (3) then determines R L . Combin- 
ing that equation with equations (4), (5), (7), and (8) to elimi- 
nate i and M x leads to an expression which can be solved 
for q numerically : 

Q 2 = R~ 2 - 17.81P 4/3 R L 2 /S /3 (l + q) 2 cos 2 d E , (9) 

where R u = RJa is a function of q and T by equation (4). 
RJa is found from R Lr and equation (3), and then i from 
equation (5). Finally, M x is determined from equation (7). 
Again, it is possible to obtain physically impossible solu- 
tions for which i is too small or for which sin i > 1.0. The 
resulting inclination distribution can also be compared to 
the optical linear polarization measurement. 

In order to optimize determination of the probability dis- 
tribution for M x , the Monte Carlo simulation generated 
normally distributed random values for 9 E with mean 28.6 
and standard deviation 2?1, and for K 0 with mean 18 km 
s l and standard deviation 3 km s'*. Uniformly distrib- 
uted random values were generated between 12 and 40 R 0 
for R a , between 0.5 and 1.0 for T, and between 0.8 and 1.0 
for Cl. 

Using 10,000 simulated parameter sets, we find M x = 
2.6t\i M q . The errors represent 1 a confidence limits. 
Only 0.3% of the sets must be discarded because sin i >1.0, 
and 3% because i < 55°. However, the resulting distribu- 
tion for i is peaked at 64° with only 13% of the probability 
between 71° and 86°, within the measured 90% confidence 
interval. Generating the same Monte Carlo simulation, but 
using a narrower range for R c (14-25 Rq) favored by stellar 
atmosphere models, yields M x — 1 -6- o s Mq. The inclina- 
tion distribution is not appreciably different from the pre- 
vious case. These estimates for M x are consistent with the 
theoretically expected range for neutron star masses. 

Because the X-ray source in this system has all the char- 
acteristics of neutron stars in systems containing a pulsar 
(except for pulsations), but is underdetermined by one mea- 
sured parameter, we shall obtain estimates for the remain- 


ing parameters by assuming values for M x within the 
theoretically expected range for neutron star masses. In 
order to avoid assumptions about i and R ( „ they are 
removed from the calculation. Afterward, results for these 
parameters can be compared to the ranges discussed above. 

The definition of f a , equation (7), is combined with equa- 
tions (3)-(5) to eliminate i and R 0 and yield an expression 
involving Q, M x , q and the measured quantities 0 E and/ c : 

Cl 2 = R £r 2 [l - (1 + q)* ,3 (fo/M x ) 213 cos 2 0 £ ] . (10) 

Equation (10) may be solved numerically for q (and hence 
M 0 ) for each set of generated parameters and an assumed 
value of M x . The inclination is then found from equation 
(7), and R a from equation (5). Finally, the expected a x sin i 
(in light-seconds) is found from 

a x sin i = 0.288 - — • (1 1) 

In 

For each assumed constant value of M x , in intervals of 
0.1 M 0 between 1.2 and 3.0 M 0 , 1000 simulated parameter 
sets were generated. Figure lOu shows the central value and 
1 <r contours for M 0 as a function of M x , and Figure 106 
shows the corresponding values of R 0 . It is clear that the 
expected probable values for M 0 and R a rise as a function 
of M x . Figure 10c shows the corresponding 1 a contours for 
the inclination. The dotted lines show the 90% confidence 
interval from variable linear polarization. Smaller values of 
i are favored by the Monte Carlo simulation. The 1 c upper 
limit from the simulation is below the linear polarization 
confidence interval for M x >2.2 M Q . 

The masses of five radio pulsars are known to within 0.1 
M 0 and are all within the range 1.3-1. 5 M 0 . Uncertainties 
in X-ray pulsar masses are typically in the range 0.5-1 M e . 
The lowest upper limit, that of Her X-l, is 1.2 Mq, and the 
highest lower limit, that of Vela X-l, is 1.55 M 0 (Lamb 
1991). Some equations of state of neutron star matter imply 
a firm mass range of 1.2-1. 6 Mq. Estimates for the compan- 
ion parameters have been derived by performing a Monte 
Carlo simulation, as outlined above, with M x drawn ran- 
domly and uniformly from the range 1.2- 1.6 Mq. With 
10,000 simulated events we obtain the following 1 c results: 
M 0 = 30^^ M 0 ; R 0 — 18 ^ 2 K 0 ; * = 66l ? degrees; and 
a = 31 1 3 2 Rq. The 99% confidence interval for a x sin i is in 
the range 48-82 lt-s. The 99% confidence interval for R 0 is 
in the range 13-26 R 0 . There are no parameter sets for 
which sin i > 1.0. Only 0.6% of the sets were discarded with 
i < 55°. The inclination in 0.1% of the sets are above the 
90% confidence interval upper boundary of 86°, and 77% 
are below the lower boundary of 71°, with 11% below 60°. 

Whether or not an input range is assumed for M x the 
most probable values of the inclination derived from our 
simulations tend to be somewhat less than the lower bound- 
ary of the measured confidence interval. Repeating the last 
simulation, using the same input parameter distributions, 
except for substituting a broad uniform distribution for 9 e 
(23°-37°), yields almost identically the results obtained pre- 
viously. This is an indication that the true inclination is 
probably close to or less than the low end of the measured 
confidence interval. More accurate measurements of the 
input parameters are needed to resolve this issue. 

Table 3 lists the values of measured parameters and 
ranges of assumed parameters used in this study. 
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Fig. 10.— System parameters as a function of X-ray mass. Dashed 
Hmitf AM° W Centra Values and S0lid curves re P resent 1 ° lower and upper 

X rav fo'rrrrr SOlar UnitS - {a> Compan,on mass a < a functio^f 
, * Com P anlon radius as a function of X-ray source 
mass, (c) Orbital inclination as a function of X-ray source mass. The dotted 
hnes enclose the 90% confidence interval from variable optical linear 
polarization measurements. 


4. ORBITAL PERIOD CHANGE 

nivl* 6 it* 1 evi dence for a change in the orbital period of 4U 
l /UO-37 came from combining the mideclipse epoch mea- 
sured by EXOSA T (Haberl et al. 1989) with those measured 

i y S l J ° neS et aL 1973 > and Copernicus (Branduardi et 
al. 1978). By separately applying the orbital profile fitting 
Sf® described earlier to the first and second half of the 
BATSE data set (about 500 days m each set), two mideclipse 
epochs have been derived from BATSE data. Fitting these 
epochs and their uncertainties together with epochs T days 

trom earlier experiments to the quadratic function 


T a =T 0 + Pn + jPP n 2 


( 12 ) 


results m a measurement of the orbital period P in days and 
period derivative P at epoch T 0 , where the independent 
variable n is an integral number of binary orbits This 
yields T 0 = JD 2,440,000 + 8900.873 + 0.015; /> = 3.41 1581 
— . 0.000027 days; and P/P = -(3.3 + 0.6) x 10~ 6 y r -* 
with x = 5.3 for 6 degrees of freedom. A linear fit to the 
same data gives X 2 = 159 for 7 degrees of freedom— an 
unacceptable fit. 
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portion of the curve implied by this result, are plotted in 
Figure 11. Here we point out that the Copernicus mideclipse 
epoch uncertainties were not explicitly quoted by the 
Copernicus experiment, and the present calculation uses 
uncertainties which give a reduced y 2 of approximately 
unity in a linear fit to their data. The epochs and uncer- 
tainties used are listed in Table 2. 

These results are consistent both with the EX OS AT 
result (P/P = -(4.0 ± 1.4) x KT* yr 1 ) and with a deter- 
m ' n / 1 t ' on th e Granat / WATCH experiment (P/P = 

(4.13 + 0.19) x 10 6 yr 1 (Sazonov et al. 1993) 
However, the uncertainty quoted for the BATSE result is 
the formal statistical error resulting from the quadratic fit 
The uncertainty quoted in the WATCH paper was deter- 
mined by repeating the fit while varying the epoch of mid- 
eclipse by up to 50 orbits. In our view, for good fits the 
statistical error resulting from the square root of the fit 
covariance-matrix element represents the minimum accept- 
able quotable error for quadratic fits such as these Thus we 
regard the error quoted by WATCH as an underestimate 
for the data plotted in Figure 11. Excluding BATSE data, 
we find a statistical uncertainty of 1.0 x 10 -6 yr" 1 for the 
WATCH observation, which has a 4.1 a significance. Inclu- 
sion of BATSE data improves the significance of the orbital 
period change to 5.5 a. 


TABLE 3 


4U 1700-37 System Parameters Used 
or Determined in This Study 1 


Parameter 

Value 

«« 

28?6 + 2?I 



18 + 3 km s’ 1 

fo 

0.0021 + 0.0010 

e 

<0.01 


r 

o.5 <r< i.o 

Q 

o.8 <n< i.o 

M x 

2-6-M M e 

M 0 

30-7* M e 



i 

661] degrees 

a 

31-1 ^0 

a x sin / b 

48 < a x sin i (lt-s) < 82 


* The uncertainties are at the 1 o level. 
b At the 99% confidence level. 
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Fig 1 1 .-Delay in mideclipse epochs with respect to a constant orbital 
period. The solid line represents the quadratic portion of the best fit to the 
epochs. 

5. DISCUSSION 

The binary 4U 1700-37 is one of the most variable per- 
sistent sources which have been observed. Using BA 15b, 
we have recorded many flaring episodes in detail and 
acquired a long-term history of the source. Two brief events 
containing clear evidence of quasi-penodic modulation 
with periods of 100 and 75 s appear in Figures 2 and Id, 
respectively. These events are similar to the observation in 
Tenma data of a 67.4 s periodic signal observed in one flare 
(Murakami et al. 1984), and the observation of four events 
containing three to four pulses each, with periods between 
96 and 248 s in EX OS AT data (Doll & Brinkman 1987). 
The signal observed by Tenma had a duration of ~ 1200 s 
and occurred at binary orbit phase 0.33, whereas the 
BATSE detections were at phases 0.16 with duration 
(100 s pulses) and 0.42 with duration at least 250 s (75 s 

PU Here we propose that these episodes may be caused by a 
beat frequency model (BFM) mechanism operating m a 
transient accretion disk near the X-ray source. In the BFM, 
the Kepler motion of a series of clumps orbiting the neutron 
star at the magnetopause accrete onto the neutron star, 
producing a modulation at the beat (QPO) frequency 
between the Keplerian rotation frequency at the magneto- 
sphere and the rotation frequency of the neutron star (Alpar 
& Shaham 1985). In this picture the magnetosphenc radius 
contracts as the mass accretion rate (and luminosity) 
increase. This could explain the different frequencies 
observed if there is a correlation between the QPO fre- 
quency and the luminosity. Because there is no permanent 
disk structure in this source, these events are only observed 
during intervals when sufficient material accumulates just 
outside the magnetopause. Although BFM was first applied 
to low-mass X-ray binaries (LMXBs), it has recently been 
tentatively identified during an outburst of the Be binary 
source A0535 + 262 (Finger et al. 1995). We have searched^ 
10 of the brightest flares in the BATSE database, and all of 
the data during TJD 8408-8511, but have not found evi- 
dence for additional QPO events. To test the BFM mecha- 
nism it will be necessary to observe enough events to 
establish a correlation of QPO frequency with .luminosity. 

The expectation for wind-fed sources such as Vela X- 
and 4U 1538-52 that P s oc P^ (Waters & van Kerkwijk 
1989) would imply a spin period for 4U 1700-37 in the 
range of P s 150-500 s. The lack of detection of pulsations 


implies either that the pulse period is very different from this 
expectation or that there is some other reason why pulses 
are not observed. This could be due to a chance alignmen 

of the magnetic and rotation axes. 

The long-term monitoring of 4U 1700-37 by the 
BATSE Earth occultation technique presented here has led 
to an improved measurement of the eclipse duration and 
ephemeris of this system. In interpreting these results we 
have used the restriction 0.8 < £1 < 1.0 derived from el lp- 
soidal brightness variations observed in optical light curves. 
This constraint has not generally been used in previous 
analyses. For example. Heap & Corcoran (1992) assumed 
i = 90° and R c = 17 R 0 . They noted that a value for M 0 ol 
27 Me is a factor of 2 lower than the masses of single stars 
of similar spectral type and position on the H-R diagram. 
The low mass implies stellar wind parameters which are 
different from those of single O stars, and therefore either 
the mass of the companion must be a factor of 2 larger or 
winds in high-mass X-ray binaries (HMXBs) af c different 
from those of single O stars. Those authors suggest that this 
problem can be resolved if the eclipse semiangle is close to 
30°. Their other assumptions lead to M 0 - 52 M Q , q > 

and a = 2 R 0 - Their scenario implies a small value ol II less 
than 0 77. Using Monte Carlo error propagation tech- 
niques, and assuming a neutron star mass consistent with 
those in other systems, we have shown that the question 
of whether the O star is undermassive cannot be hrmly 
resolved. However, the data do favor smaller masses. 

An orbital period change has been observed for only two 
other HMXBs: Cen X-3 (Nagase et al. 1992) and SMCX-1 
(Levine et al. 1993). The system parameters and orbital 
period change determined in this paper can be used to 
assess the cause of the orbital period change. This change 
could be caused by wind-driven angular momentum loss, or 
bv tidal interactions, or by a combination of both. 

Our results can be compared to a recent ballistic 
(numerical) calculation of the contribution of wind-driven 
angular momentum loss to orbital period changes in high- 
mass binary systems. The relation between the rate of 
orbital period change, the mass-loss rate of the companion, 
the average angular momentum loss per unit mass of win 
particles in the center of mass of the binary system (h cm ), an 
the mass ratio q is (Brookshaw & Tavani 1993, hereafter 

BT): 


M- 


p 

p 


< _ i 

cm 3 1+ q 


Mo 

M 0 


Tl ° , (13) 


where the mass fraction of the wind captured by the 
primary is assumed to be negligible for this purpose 

Tables of h cm as a function of q, £2, and wind velocity V 
normalized to the orbital velocity of the secondary around 
the center of mass, V cm = aa>/(\ + q), have been compiled by 
BT. From our results we may take F cm - 20 ± 3 fcm s • 
The escape velocity of the wind is 617 km s , and the 
terminal velocity (Heap & Corcoran 1992) is 2100 km s , 
which imply normalized velocities V - 31 + 5 ana 
y — j06 + 28, respectively. From Table 1 of BT, at large 
velocities (F > 30) h cm apparently approaches a nearly con- 
stant value, and changes negligibly with velocity but 
increases with q. From Table 13 of BT, it ,s apparent that 
h m also increases with £2. For q = 30, 0=1. BT find C m 
9.5. If we also use M a = -6 x 10 6 M 0 yr 1 Haberl et al 
1989) and M a = 42 A# 0 , then equation (13) yields an orbital 




period change P/P = -3.8 x 10 6 yr" \ within 1 a of our 
measurement. Thus we cannot rule out the possibility that 
the period change is due entirely to angular momentum loss 
in the stellar wind. However, we note that the results pre- 
sented in BT assume corotation of the companion with the 
binary orbit. Recently, h cm has been calculated as a function 
of T (M. Tavani 1995, private communication). In general, 
h cm increases with T. The smallest possible value of h cm for 
this system, at q = 15, Q = 0.8, T = 0.5, is h cm = 1.5. 
Assuming also M 0 — 20 M 0 implies an orbital period 
decrease due to mass loss which is 25% of that observed. 
Allowing also for a ~50% uncertainty in M 0 , we set a 
lower limit of 10% on the contribution of mass loss to the 
observed orbital decrease. Thus, considering also the 
q = 30, Q — 1, T = 1 case discussed above, we conclude 
that ballistic calculations imply that between one-tenth and 
all of the observed orbital period decrease is due to mass 
loss by the companion star. 

As an additional complication we note that these ballistic 
calculations do not include the radiative forces acting on 
the wind. How this impacts the calculation of wind-driven 
orbital angular momentum loss depends on the details of 
the accretion. The velocity law which is usually assumed is 
of the form v = 1^(1 — R 0 / r f * where p, x , is the terminal wind 
velocity, a is the wind acceleration parameter, and r is the 
distance from the center of the secondary star. This wind is, 
however, disturbed by the presence of the X-ray source, 
which photoionizes surrounding material, decreasing the 
radiative cross section and stopping the radiative acceler- 
ation in some region (Stromgren sphere) around the X-ray 
source. Haberl et al. (1989) fitted EXOSAT data to the 
observed column density as a function of phase. They found 
a = 0.43 and a dimensionless ionization parameter £ = 190, 
which measures the size of the Strdmgren sphere. With 
these values, the wind reaches escape velocity at about r = 
1.1 R 0 , while the Stromgren sphere does not begin until 
about r= \AR 0 . Thus, if these parameters are approx- 
imately correct, neither the radiative acceleration nor the 
X-ray photoionization should have an appreciable impact 
on the mass-loss calculation. 

The contribution of tidal interactions to the orbital 
period change can be estimated by computing the circular- 
ization time for a radiatively damped tide (Lecar, Wheeler, 
& McKee 1976, eq. [19]). For the central measured param- 
eter values of 4U 1700 — 37, we find that t circ = 7 x 10 6 yr. 
Propagating 1 a errors on the measured quantities shows 
that the tidal mechanism can supply 4.3 tfg percent of the 
orbital period change. That is, the 1 o upper limit is about 
12.4% and the central value is about 4.3%. The 95% con- 
fidence upper limit is 21%. This is consistent with the above 


result that the period change due to mass loss is between 
10% and 100% of the total. Thus it appears likely that in 
this system mass loss is the dominant cause of the orbital 
period change. 

6. CONCLUSION 

Summing almost 1000 days of BATSE occultation data 
has permitted determination of the eclipse semiangle of 4U 
1700—37 with a 1 a uncertainty of about 2°. This measure- 
ment, combined with the optical and UV measurements of 
the radial velocity semiamplitude, a range constraint on the 
Roche lobe filling factor, and the assumption of a circular 
orbit have led to a new determination of the system param- 
eters. Propagation of errors on experimental quantities 
shows that the mass of the X-ray source is very probably 
within the theoretically expected range for neutron stars. 
The question of whether the companion is undermassive 
cannot at this time be firmly resolved, but the evidence does 
favor smaller masses. No evidence for a coherent pulse 
period between 2 and 700 s has been observed. 

An OTTB model spectrum with a temperature of 25 keV 
is consistent with BATSE data, independent of source 
intensity or orbital phase. Observations of flaring behavior 
show variations on timescales from 30 s to more than 1 day. 
The distribution of flare sizes is consistent with an exponen- 
tial model. Transient events of a quasi-periodic nature 
observed in BATSE data as well as by other experiments 
hint at the possible formation of transient disk structures. 

Combining two mideclipse epochs determined from 
BATSE data with those determined from earlier experi- 
ments leads to a value for the rate of orbital period decrease 
which is significant at the 5.5 a level. The contribution of 
mass loss in the companion wind to this decrease is between 
10% and 100%, while the contribution from tidal inter- 
actions is probably less than 20%. 

BATSE occultation observations of 4U 1700 — 37 over 
the next 5 years will permit a further reduction of the uncer- 
tainty in 0 E and continued monitoring of the orbital period 
change. Accurate measurements of 0 E by other experiments 
and improved measurements of other parameters are 
needed to sharpen the determination of the system param- 
eters. In addition, higher time resolution measurements of 
the orbital cycle could do more to reduce uncertainty in the 
orbital period change. 

We thank Jan van Paradijs for useful discussions about 
the system parameters, and Marco Tavani for valuable dis- 
cussions about the interpretation of the orbital period 
change. 
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